Probing the 3.6 Micron CIRB with Spitzer in 3 DIRBE Dark Spots 



L. R. Levenson & E. L. Wright 

Department of Physics and Astronomy, University of California, Los Angeles, CA 

90095-1562 

levensonOphysics . ucla . edu 
ABSTRACT 

We observed three regions of the sky with Spitzer in which the Cosmic In- 
fraRed Background (CIRB) has been determined at 3.5 /im using the method of 
subtracting 2MASS stellar fluxes from zodiacal light subtracted DIRBE maps. 
For each of these regions we have obtained 270 seconds of integration time per 
pixel with IRAC on Spitzer over the central square degree. We present galaxy 
counts in each of these approximately 1 square degree IRAC surveys. Along 
with deep galaxy counts in the Extended Groth Strip and GOODS North, we 
are able to compare the galactic contribution to the CIRB with the "DIRBE 
minus 2MASS" determined L-band CIRB. Using the profile-fit photometry pack- 
age GIM2D, we find a substantially larger flux contribution to the CIRB than 
that determined using aperture photometry. We have also made the first rigor- 
ous analysis of the uncertainties in determining the CIRB via galaxy counts in 
Spitzer images using a Monte Carlo Markov Chain simulation of our data anal- 
ysis. Using a simple broken power law model for galaxy counts as a function 
of magnitude we find a most probable contribution to the CIRB from galaxies 
at 3.6 /xm of 10.81?;} kJy sr-i (D.O^J;^ nW m ^ sr Even with this restricted 
model, however, we find that galaxy counting does not strongly constrain the 
CIRB from above. We are able to find solutions in which the CIRB runs away 
to large intensities without the need for an additional diffuse source. 

Subject headings: cosmology: observations — diffuse radiation — infrared: galax- 
ies 



1. Introduction 



The Cosmic InfraRed Background (CIRB) is the aggregate of the short wavelength 
radiation from the era of structure formation following the decoupling of matter and radiation 
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in the early universe, redshifted to near infrared wavelengths by cosmological redshifting and 
far infrared wavelengths by dust-reprocessing, i.e. absorption and re-emission of starlight 
by intervening dust. Particle decay models also allow for a contribution to the CIRB from 
photons p roduced in the de cay of weakly interacting massive particles such as big bang relic 
neutrinos (IBond et al.lll986l ). Hence, the cumulative light from galaxies is a strict lower limit 
on the CIRB. A determination of the total contri bution of re s olved galaxies to the CIRB 
was made recently via galaxy number counts by iFazio et al.l (120041. Using the InfraRed 
Array Camera (IRAC) on the Spitzer Space Telescope (lEisenhardt et al.l |2004| ). surveys 
were done of the Bootes region, the Extended Groth Strip and a deep image surrounding 
the QSO HS 1700-1-6416. After integrating the light from galaxies from the 10th to the 
21st magnitudes, a total integrated intensity of 6.5 kJy sr~^ (5.4 nW m~^ sr~^) at 3.6 
nm is reported. This i s less than ha l f of t h e CIRB determined at t hat wavelength by 
Dwek fc Are^ JlQQsh. ICoriian et all J200oh . IWright fc Rees3 J200oh . IWright fc Johnson 



(j200l|) and 



Levenson et al 



J2007h which directly measured the CIRB via the subtraction of 



foregrounds from all sky maps from the Diffuse Infr aRed Background Ex per iment (DIRBE) 



on th e COsmic Background Explorer (COBE, see iBoggess et al.l (119921 )). iLevenson et al. 
J2007h determine a 3.5 fim CIRB of 15.6 ± 3.3 kJy sr^^ (13.3 ± 2.8 nW m'^ sr^^). This leaves 



a greater than 2a gap between the direct measurements of the CIRB and the lower limit from 
galaxy counts. To investigate this discrepancy, we have obtained Spitzer observations of the 
central square degree of three of the "dark spots" used in the "DIRBE Minus 2MASS" direct 
determi i iations of the CIRB, where " DIRB E Mi nus 2MASS" ref e rs to the series of papers by 
Wrightl J2001ah . IWright fc JohnsonI J200lh and Levenson et al.l (120071 ) which estimated the 
CIRB intensity by subtracting the IWrightl (119981 ) zodiacal light model, 2MASS stellar fluxes 



and a faint star count mod el from the t otal D IRBE intensit y. The galaxy flux cont ribution 
to the CIRB determined by lFazio et al.l (120041 ) uses Bertin's (IBertin fc Arnoutslll996l ) Source 
Extractor (SE) for both source detection and aperture photometry. While we also use SE for 
source detection, we employ a profile-fit photometry method using the IRAF package GIM2D 
(ISimard et al.ll2002l ) in an attempt to include more of the faint fringes of faint galaxies in 
our cumulative count of the light from these galaxies. 

To this point, only lower limits from galaxy counts have been reported in the literature. 
Here we have estimated uncertainties resulting from counting (Poisson) and cosmic variance, 
as well as from the flux measurements, whether from profile-fit or aperture photometry. At 
faint magnitudes (m > 16) the sky becomes so crowded that confusion becomes a significant 
barrier to flux estimation using either method. Errors become more pronounced at faint 
magnitudes where confusion is significant, but even at bright magnitudes profile fitting tends 
to overestimate the flux, while aperture photometry underestimates it. We have corrected 
for these errors using a rigorous analysis employing a Monte Carlo Markov Chain (MCMC) 
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simulation of our analysis by constructing a simple broken power law model for the "real" 
galaxy counts and using the Markov Chain to explore the parameter space of models that 
could lead to our data. Two simple priors were used, one requiring that the integrated flux 
converge and another weighting the bright end galaxy counts near Euclidean counts. We are 
now able, drawing from the posterior distribution, to compute the CIRB intensity with the 
first robust error bars on the galactic contribution to the CIRB. 



2. Data 



Our most recent mea surement of the L-band CIRB (ILevenson et al.l 120071 ) used the 
method of IWrightl (j2001al ) to subtract stellar fluxes from the 2 Micron All Sky Survey 
(2MASS) Point Source Catalog (PSC) from zodiacal light subtracted DIRBE maps in 40 
regions of the sky. All 40 of these regions are at high galactic latitude and scattered widely 
in ecliptic latitude. The three regions of the sky chosen for observation with Spitzer were 
selected from these 40 DIRBE "dark spots." The regions chosen are centered at widely dif- 
ferent ecliptic latitudes of +20 (region 2), +50 (region 1) and +70 (region 3) degrees. Actual 
region coordinates are shown in Table [H The observing strategy was to use the mapping 
capabilities of IRAC to map the central square degree of each of these regions with a total in- 
tegration time of 270 seconds per pixel. Observations were divided into 3 sets of 3-30 second 
frames at each pointing, always at a new point in the dither table. Thus, each pointing was 
visited on nine separate occasions with three observations approximately 45 seconds apart 
as the dithered map was executed within each AOR, and returning three times to repeat 
this process anywhere from hours to days to months later depending on the scheduling of 
separate AOR's for observation. This redundancy allowed robust removal of artifacts such 
as cosmic rays, whie the dithering ensured that a bad pixel would be backed up with several 
observations with good pixels. Once all of the individual observations were mosaiced, the 
total usable coverage obtained in each region is 0.60, 0.84 and 0.82 square degrees in regions 
1, 2, and 3 respectively, where region 1 is smaller due to a long gap between observations 
that allowed a rotation of the telescope. 

As shown in Figure [H source counts in these three regions begin to be severely incom- 
plete between the 16th and 17th magnitudes. Galaxies at these magnitudes are the peak 
contributors to the CIRB. In order to obtain counts at magnitudes deeper than the peak 
contribution to the CIRB, the use of deeper images was required. Counts at intermediate 
depths were done using a deep image of the Extended Groth Strip (EGS) observed as part of 
the GTO project "The IRAC Deep Survey." The deepest counts were done using the Great 
Observatories Origins Deep Survey (GOODS) North field (IDickinson et al.ll2002l ). 
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3. Analysis 



Basic Calibrated Data (BCD) images were combined into mosaics in each region using 
MOPEX. Source detection in each regio n, as well as in the dee p image mosaics, was then 
performed using Source Extractor (SE) ( Bertin fc Arnout^ 1996 ). A conservative detection 
limit required 5 adjacent pixels at least 5a above the SE determined background. Aperture 
magnitudes for ea ch source w ere also determined using SE's mag_auto which uses an ad- 
justable elliptical (IKronl 119801 ) aperture. However, since aperture photometry can miss the 
faint fringes of galaxies, especially in crowded fields where large apertures can not be used, 
a more s ophisticated metho d of profile fit photometry was employed. The IRAF package 
GIM2D (ISimard et al.ll2002l ) fits a 12-parameter, bulge plus disk profile to each source and 
determines the photometry from the resulting model. A comparison of the aperture and 
profile-fit galaxy-only counts in each region are shown in Figure [H Removal of stars from 
the total counts is described in detail below. The resulting intensity contribution of these 
galaxies to the CIRB using the two methods are shown in Figure [21 Profile fitting results 
in magnitudes approximately 0.5 mag brighter, on average, for the same source. This re- 
sults in approximately 1.6 times more intensity contributed to the CIRB. All magnitudes 
quo ted here are Vega -relative magnitudes where Vega has a flux density at 3.6 /im of 277.5 
Jy jPazio et al.ll2004h . 



In order to determine the contribution of galaxies to the CIRB, stars must be removed 
from the raw source catalogs. At brighter magnitudes, stars and galaxies can be separated 
morphologically. GIM2D returns a half-light radius for each fitted source. This half-light 
radius is as reliable an indicator of a galaxy as visual inspection of the galaxy image, for 
galaxies brighter than magnitude 13.5. At magnitudes fainter than 13.5, both visual inspec- 
tion of the source image and the GIM2D half-light radius become ambiguous indicators of 
morphology and another method of star/galaxy separation is required. Regions 1 & 2, as 
well as th e EGS and GOODS-N are covered by the Sloan Digital Sky Survey (SDSS) data 
release 5 (lAdelman-McCarthy. J. et al.l 120071 ). All of the detected sources in the resulting 
SE/GIM2D catalogs from these images were fed into the SDSS/Skyserver Imaging CrossID 
web tool. Approximately half of the sources in our IRAC images had matches within 3" in 
the SDSS catalog. Any source (m > 13.5) with a match within 3" classified by SDSS as a 
star was removed from the catalog. All sources classified by SDSS as galaxies, along with 
all sources red enough to be detected in our IRAC images and not detected by SDSS were 
considered to be galaxies and included in the total intensity of the CIRB. The nearby 8th 
magnitude galaxy NGC 5927 in region 1 was not included in the total CIRB estimate. 

As region 3 was not covered by SDSS, it was necessary to remove stars with m > 13.5 
statistically using the average of the star/galaxy ratios at each magnitude in regions 1 and 
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2 multiplied by the ratio of the sines of the galactic latitudes. This assumes a csc|b| scaling 
for star counts. These ratios were applied to the raw source counts at each magnitude in 
region 3 using, 



\ _Q^' f 'f^staA / sin&2 \ ^ / ristaA f sin hi \ 



(1) 



where hi is the galactic latitude of region i. Resulting galaxy-only counts in all three regions 
are shown in Figure [3l Hereafter, all counts refer to those determined using profile-fit mag- 
nitudes. While the fractional variance 5N / (N) in the galaxy counts per unit solid angle at 
a particular magnitude is as high as 11%, 6N/{N) of the total counts per unit solid angle 
between the three regions is less than 1%. The fact that the total galaxy counts do not show 
a dependence on galactic latitude, b, indicates successful removal of stellar sources in these 
three regions and supports the above assumption of a csc|b| scaling for the star counts. 

Figures [3] and m clearly show that counts in regions 1-3 are incomplete by m = 16.5. 
Since this is near the peak contribution to the CIRB, accounting for this incompleteness is 
crucial to obtaining an estimate of the galactic contribution to the CIRB. We chose to use 
deeper images to obtain real counts at these important magnitudes rather than rely on an 
estimate of the incompleteness in our surveys. Figure H] shows the average of the counts in 
our three regions along with the deeper counts, obtained using the exact same methods, from 
the EGS and GOODS-N. Final observed counts are then the averaged counts from our three 
regions for m < 16.0, averaged counts from the EGS and GOODS-North for 16.5 < m < 
17.5, and GOODS-North for m > 18.0. 

Completeness was estimated using the standard Monte Carlo method of inserting sample 
galaxies of known magnitudes into the science images at known coordinates and re-running 
the source detection. The sample galaxy images used were model images created by GIM2D 
in the process of fitting actual galaxies. Sources were considered recovered if a new source was 
detected within ±1.5 pixels of the insertion coordinates, or, if an old source was brightened 
by at least half of the inserted galaxies flux. In the shallow fields, for magnitudes less than 
16.5, at each half magnitude, 200 copies of at least 31 different model galaxies were inserted 
and re-detection attempted in 4 iterations making a minimum of 24800 sample sources 
inserted. At fainter magnitudes, since the deep fields are smaller and more crowded, only 40 
copies of each galaxy were inserted in a single iteration. However, at least 34 different model 
galaxies were used and 10 iterations of the process were performed resulting in a minimum 
of 13600 sample sources inserted at each half magnitude down to m = 19.5. Completeness 
is already below 40% at m = 19.5, so observed counts were cut off at this magnitude. 

Completeness estimates for magnitudes 10.5 - 19.5 are shown in Table El Our an alysis 



finds the deep images to be substantially less complete than the iFazio et al.l (120041 ) esti- 
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mates, even in the same image (EGS). Several effects may be causing this difference. The 
first is the detect ion hmit used. While SE was used for source detection in both cases, 



Fazio et al.l (120041 ) use a 2.5a detection limit while we have used a more conservative 5a 
limit. Thus, given the same spatial scale, they will detect a galaxy with a surface brightness 
twice as faint, and thus go 0.75 magnitudes deeper. In addition, the magnitude at which the 
completeness is determined is given here by our profile-fit magnitudes as opposed to their 
aperture magnitudes. We have shown here that our magnitudes are, on average, approxi- 
mately 0.5 magnitude brighter for the s ame source t han th ose determined with SE aperture 



photometry. Deep field photometry in iFazio et al.l (120041 ) was done with custom software 



however, so it is difficult to know just how these compare. 

The remaining difference may be due to the observational characteristics of the partic- 
ular galaxies used for completeness estimation. Here, the galaxies inserted into the science 
images were chosen at random from each half magnitude bin. This results in a wide variety 
of observational characteristics such as the peak of the surface brightness distribution, the 
ellipticity, bulge fraction, and spatial extent, all of which affect the rate at which galaxies 
of the same magnitude are recovered. For example. Figure [5] shows a sample of two of the 
galaxies used, labeled A and B, which have very different peak surface brightnesses. The 
indicated galaxies in both the right and left panels were inserted at the same coordinates 
and both have profile-fit magnitudes of m = 16.5. However, the galaxy inserted into the left 
frame (galaxy A) has a peak surface brightness of 207.0 kJy sr^^ while the galaxy inserted 
into the right hand frame (galaxy B) peaks much higher, at 259.2 kJy sr^^. And while 
galaxy A actually has a larger total flux, it is detected at only 24 of 40 inserted positions in a 
single iteration compared with 35 of the 40 insertions for galaxy B. However, the peak of the 
surface brightness distribution is not the only feature affecting completeness. Figure [6] shows 
that while there is a trend with the peak surface brightness, the completeness at a given 
peak surface brightness is also widely scattered. Figure [7] shows that the way the remaining 
flux is distributed about the peak, affected by the galaxy's inclination and morphology, also 
plays a role. The left hand frame is the same as that in Figure O where galaxy A has a 
peak surface brightness of 207.0 kJy sr~^ and a 60% detection rate. The galaxy inserted into 
the right hand frame (galaxy C) has a similar peak surface brightness of 211.9 kJy sr~^, but 
a detection rate of 78%. Visually, it is obvious that galaxy C has a higher ellipticity than 
galaxy A (0.7 vs. 0.4, based on the GIM2D fit). Less obvious to the eye, however, is that 
the latter has a half-light radius approximately 40% larger. For this work it is not necessary 
to conduct an exhaustive investigation into the detailed variations of galaxy properties such 
as ellipticity, inclination, bulge fraction, etc. and their effect on the ability of SE to detect 
a given galaxy. What is important, however, is to note that the difference in detection rates 
of model galaxies based on a wide variety of real galaxies introduces a large scatter into the 
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completeness estimates. Inclusion of this scatter will be crucial in determining uncertainties 
in the completeness corrected number counts and will be discussed in the following section. 

3.1. Uncertainties 

The total intensity of the 3.6 fim CIRB is computed by summing over the contributions 
from galaxies binned in 0.5 magnitude bins centered at each half magnitude: 

/c//?s(3.6/im) = ^ n{m)f{m) , (2) 

m 

where n(m) is the number of galaxies per steradian at magnitude m and f(m) is the flux of 
the center of the magnitude bin. A robust estimate of the error associated with estimating 
the CIRB by counting galaxies in IRAC images using Source Extractor and GIM2D requires 
an analysis of cosmic variance in addition to the Poisson shot noise in the counts, as well 
as the error in determining the completeness of the source detection and in determining 
the flux of each galaxy. Uncertainties in the counts, completeness and flux are discussed in 
Sections 13.1.11 - Section [3. 1.31 

3.1.1. n(m) 

Uncertainties in n(m) are a combination of Poisson errors and cosmic variance. For 
m < 16.0, where n(m) is measured in our three regions, the measured fractional variance at 
each half magnitude: 

5n _ ^/{{n-{n)Y) ^ 
(n) {n) 

where (n) is the average n(m) in the three regions, should include both the Poission noise 
and the cosmic variance since these fields are widely separated, allowing us to assume there 
is no spatial correlation between the counts in each field. These measured values, however, 
in some bins fall below the Poisson error based on the counts. The values in Table 1 for 
m < 16.0, which are taken to be the uncertainties in those counts, are then the larger of the 
Poisson error or the measured fractional variance. 

In the deep fields (EGS, GOODSN), at fainter magnitudes, cosmic variance is compa- 
rable to the Poisson noise. While the fields are small in angular size, they are deep pencil 
beam surveys and thus sample a large volume of space. We can use a geometrical argument 
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to estimate the cosmic variance in the deep fields. The average of the SDSS photometric 
redshift estimates for galaxies with matches in the SDSS catalog is z ~ 0.5 These images 
probe much deeper than z = 0.5, however, we can use this to put an upper limit on the 
cosmic variance. Taking the average depth z ^ 0.5 and a ACDM cosmology with = 0.27, 
Q\ = 0.73, we reach a proper radial distance of ~ 2700/i~^ Mpc {h = Ho/100 km s~^ Mpc~^). 
Breaking this up into pieces of length 8h~^ Mpc, this spans over three hundred 8h~^ Mpc 
slices. Taking erg ~ 1, where erg is the relative density contrast in a spherical region of radius 



8h^^ Mpc, and the variance proportional to A^jje/^s where Nfi^ids is the number slices of 
length 8h~^ Mpc, this limits the fractional cosmic variance in the counts to less than about 
6%. This is consistent with the measured 6N/ (N) between the EGS and GOODS North 
fields of 5.1, 4.1 and 6.4% at m = 16.5, 17.0 and 17.5. The error due to cosmic variance in 
the deep fields is then taken to be ± 5%. The total uncertainties in the number counts per 
magnitude at faint magnitudes, used in computing in the next section, are then, 



where n is the number per magnitude per square degree, is the solid angle of the region 
in square degrees and the factor 0.5 comes from using half magnitude bins. These are the 
values shown in Table [2] for m > 18, and the quadrature sum of these uncertainties in the 



As discussed above, completeness at each magnitude was determined by inserting ap- 
proximately 40 different sample galaxy images into the science frames at 200 or 40 different 
random positions, depending on the image, and re-running the source detection and cross- 
correlating the new catalog with the old. Thus, for each sample galaxy, i, we obtain a 
completeness estimate, compi, from the ratio of the number detected to the number in- 
serted. The value of compt varies widely with i due to the different surface brightnesses and 
morphology of each sample galaxy. We then run Nj^r = 4 or 10 iterations of this process, 
depending on whether we've inserted 200 or 40 copies of our sample galaxy. For the j*'' iter- 
ation, the completeness, compj, is taken to be the mean of the estimates for each galaxy and 
the completeness uncertainty, acompj, is conservatively taken to be the standard deviation 
of the distribution of the values compi. The completeness at each half magnitude is then : 




(4) 



EGS and GOODS North is shown for 16.5 < m < 17.5. 



3.1.2. Completeness 



comp{m) = (compj) , 



(5) 
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and 



^comp yfTT') 



\ 



loiter 



iter 



a 



comp,j 



j=0 



The completeness corrected, observed number counts them become: 

nobs,corr{m) = Hobsim) / comp{m) , 

with a total uncertainty of, 



(6) 



(7) 



(m) 



2 / \ 2 

^n.obs \ ^ I ^comp \ 

^obs 



\comp J 



f nobs Y 
\ comp J 



3.1.3. f(m) 

Once the sample galaxies were inserted into the science frame and re-detected for the 
incompleteness estimates, for the last run, the entire analysis pipeline was re-run on those 
sample galaxies that were successfully detected as new sources within 1.5 pixels of the inser- 
tion coordinates. After the sample galaxies were re- fit, the new fluxes determined by GIM2D 
were compared with the actual flux of the inserted galaxy. Histograms of the total flux from 
the re-fit, with the total of the inserted pixel values subtracted, are shown for each half 
magnitude from m = 14 to m = 18 in Figure [HI The histograms are normalized to a total of 
one. Vertical lines in the plots show the widths of the 0.5 magnitude bins. One can see from 
the histograms that the modes of the distributions are near zero, i.e. the most likely value of 
the re-fit flux is the actual flux. However, we also see that many of the new fluxes spill out of 
the original bins into neighboring magnitude bins with a clear skew toward overestimation 
of the flux. Given this information, however, we use a Monte Carlo Markov Chain (MCMC) 
simulation of our analysis to put constraints on the possible real distributions of n(m) that 
can give rise to counts we observe. 

The model for the counts of galaxies as a function of galaxy flux, n(f) [mag~^ deg~^], is 
a simple broken power la^^|: 

= (1 + (Mfr) • 

where A, a, (3 and /o are adjustable parameters and the bright end counts are Euclidian if 
a = 1.5. Two prior assumptions were made. The first is that the integrated flux converges. 



^In principle, one can always generate a very large CIRB by adding terms to n{f) that give large counts 
for unobservably small fluxes. Our simple form for n{f) rules out this possibility a priori. 
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achieved by requiring a faint end slope a — (3 < 1. The second is a prior on the bright end 
slope, a, such that before subjecting the model to the data, we assume a = 1.5 ± 0.25. The 
likelihood function is given by L = e~*^'^-^^ where 



(10) 



The posterior probability is then given by: 

P{A, a, /3, fo\nobs,corrim)) = e"°-^^V(« - l3)p{a) 



where 



and 



p{a - (3) 



1, for a — /3 < 1 
0, for a - /? > 1 



p{a) = exp —0.5 



a — 1.5 
0.25 



(11) 
(12) 

(13) 



The p arameter space is then explored using Markov Chains (see ISmithI (1199 ih or [Binder 
(1l997l )). starting with a various sets of initial parameters. The MCMC simulation proceeds 
as follows: 



1. Compute the modeled "real" counts at each half magnitude for 10.5 < m < 19.5, given 
the current set of parameters using Equation [9l 

2. Compute the modeled "observed" counts by multiplying the modeled "real" counts 
vector njnodoi.cai t>y the transfer matrix, h, and summing the resulting counts at each 
flux, j, i.e. n„iodei,ts,j = 'E^i^modeUsAij- The i*^ row of the matrix, one for each half 
magnitude bin, is given by the histograms in Figure [8] with the actual flux at the bin 
center added back in. Histograms for m < 14 and m > 18 are the m = 14 and m = 18 
histograms respectively, scaled by the central flux of the appropriate bin. 

3. The resulting modeled "observed" counts are spread over the full range of fluxes, j, 
and must then be binned into the same 19 half magnitude bins as the observed counts. 
We can then compute the posterior probability as given by Equation [TTl 

4. Generate a random change in the parameters 6 A, 6a, 5(3 and 5fo- 

5. Repeat Steps 1-3 

6. If the probability of the new parameters, Pnew, is higher than that of the old parameters, 
Poidi the new parameters are accepted and the process is repeated from Step 4 
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7. If 

Pnew ^ Poidi the new parameters are accepted with the probabihty Pnewj Poid snd the 
process is repeated from Step 4. This facihtates exploration of the entire parameter 
space by sometimes accepting a downhill step. Otherwise, the new parameters are 
rejected and the process is repeated from Step 4. 

Several chains were run with 200,000 steps each. 

4. Results 

Figure [9] shows the posterior probability distribution and the resulting CIRB for the 
models with the highest probability in each bin. The distribution gives a most probable 
contribution from galaxies to the CIRB of 10. 8t^ } kJy sr~^ (O.Olg 9 nW m~^ sr~^) where the 
asymmetric range is given by the full width of the histogram where the probability falls from 
its maximum by a factor e~3. Parameters for the most probable model are given in Table [31 
One can see that even given the prior on a, the most probable bright end slope for the 
counts is steeper at a = 1.85. Figure [TO] shows the 68%, 90% and 99% confidence limits on 
the parameters a and f3. The dashed line shows the convergence prior. Figure [TT] shows the 
same confidence limits on the CIRB intensity at m = 16.5 and the integrated CIRB. We can 
see that the intensity of galaxies around m = 16.5 is well constrained, while the integrated 
CIRB intensity is not on the bright end. Figure [12] shows the same confidence limits on the 
faint end slope in the power law model of the counts {a — (3) and the magnitude at which the 
break occurs (Break Magnitude = — 2.5 log(/o/-Fo)) where = 277.5 Jy is the fiux density 
of Vega. Indicated by the cross in these three figures are the most probable values for 
those parameters. Counts and intensities as a function of magnitude derived from the most 
probable set of parameters are shown in Figures [13] and [H] Figure [13] shows the completeness 
corrected observed counts, the model "observed" counts and the model "real" counts, versus 
magnitude. The CIRB intensities resulting from the raw galaxy counts, the completeness 
corrected counts and the most probable model are shown in Table [U Figure [H] shows the 
modeled "real" counts and the resulting contribution to the CIRB at each magnitude. 

4.1. Aperture Photometry 

As a check, we ran the same MCMC analysis using counts obtained with SE aperture 
photometry. Aperture fiuxes of the galaxies inserted at each magnitude for the completeness 
testing were compared to the total of the pixel values actually inserted for each source and 
binned as before. The resulting comparison is shown in Figure [15] where the magnitude 
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at which the histogram is centered is shifted by one half-magnitude with respect to those 
in Figure [H This shift, along with a similar shift in the completeness estimates, is due to 
the use (in measuring the completeness and flux uncertainties) of galaxies selected by their 
profile-fit magnitudes. We have seen that these are, on average, 0.5 magnitude brighter 
than the SE aperture magnitudes. We can see from Figure [15] that aperture photometry 
tends to underestimate the flux. The MCMC analysis was run again with these histograms 
used in place of those in Figure [8] to scatter the modeled "real" counts. The resulting 
probability distribution for the intensity of the CIRB is shown in Figure dni Here we have a 
somewhat tighter distribution, but we still have a long non-Gaussian tail extending to larger 
intensities. The full width of the histogram at P = e'^Pmax gives a CIRB of 9.2^07 kJy 
sr~^ nW m~^ sr~^) and the best fit parameters are shown in the second column of 

Table [31 The resulting model "observed" counts, along with the model "real" counts and the 
completeness corrected observed aperture photometry counts are shown in Figure [T71 The 
CIRB intensities resulting from the raw aperture photometry galaxy counts, the completeness 
corrected aperture counts and the resulting most probable model are shown in Table [H The 
remaining la difference between the aperture and profile-fit intensities is approximately what 
one would expect, given that the two measurements are not independent, and is due mainly 
to the difference in the faint end slope of the model counts. This slope is determined by the 
most uncertain faint counts as seen in the middle panel of Figure [181 However, as shown in 
the bottom panel of Figure [181 the MCMC analysis does give a more consistent picture of 
the counts determined by the two methods. We report here the profile-fit results, over those 
from aperture photometry, based on the ratio of the maximum probabilities in each case: 
[aper) / Pmaxiprofile) fa 0.04. 



Discussion 



Wrightl (l2001bl ) extrapolated the SDSS determined optical luminosity function into 



the NIR and fo u nd it to be a factor o f 2.3 h igher than that determined with 2MASS by 
Kochanek et al.l (l200ll ) and I Cole et al.l (l200ll ). Similar discrepancies of factors of 1.5 and 
2.1 were also seen among optical luminosity functions when comparing the SDS S luminosity 



funct ion with those determined using the 2 deg ree Field Galaxy Redshift Survey (iFolkes et al. 



200ll ) and the Las Campanas Redshift Survey flLin et al.lll996l ). iBlanton et al.l (120011 ) showed 
that the optical discrepancy was explai ned by the use o f larger apertures by the SDSS. The 
NIR discrepancy was not explained, but IWrightl (l2001bl ) suggested that inclusion of the faint 
fringes of galaxies using larger apertures in the NIR could be at least a partial explanation for 
the discrepancy seen between the SDSS extrapolated and 2MASS luminosity functions. The 
Wrightl (l2001al ) DIRBE minus 2MASS subtraction, recently extended to 40 new regions of 
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the sky in lLevenson et al.l (120071 ) resulted in a detection of a statistically significant isotropi c 



background of 15.6 ± 3.3 kJy sr'^ (13.3 ± 2.8 nW m'^ sr'^) at 3.5 /im. iFazio et al.l (120041 ) 



using custom aperture photometry, determined a total contribution of resolved galaxies to 
the CIRB of only 6.5 kJy sr~^ (5.4 nW m~^ sr~^), less than half of that value. Hence, we 
set out to include the the faint fuzzy fringes of the galaxies. Aperture photometry is limited, 
especially in crowded fields where the edges of the galaxies overlap and l arge apertures can 
not be used. Here we have employed profile- fit photometry and raised the lFazio et al.l (120041 ) 
hard lower limit and replaced it with a Gaussian lower limit using counts, obtained with 
GIM2D, of galaxies in 3 DIRBE Dark Spots, as well as deep images of the Extended Groth 
Strip and the GOODS North field. We are able, using the machinery of MCMC simulation, 
to obtain a model for the distribution of galaxies, n{m), that gives a most probable total 
galactic contribution to the CIRB of lO.StH kJy sr'^ (Q-Otli nW m ^ sr ^). These values 
are consistent, for the first tim e, with prev i ous d eterminations of the CIRB in the L band 
shown in Tabled excepting the lFazio et al.l (120041 ) lower limit, which, based on the la lower 
bound from the width of the probability histogram at P = Pmaxe"^'^, is 3.5cr below the most 
probable intensity determined here. 

A third constraint on the CIRB comes from th e attenuation of Te V 7-rays by e+e^ pair 
production through collisions with CIRB photons. iDwek et al.l (|2005l ) derived the attenuat- 
ing effect on 7-r ay photons f r om th e Blazar PKS 2155-304 of the l/xm spike in the CIRB 
suggested by the iMatsumotol ( 120051 ) IRTS measurements shown in Figure [191 That analysis 
showed that such a spike is unlikely to be of extragalactic origin and that the CIRB can not 
be significantly higher than the intensity deriv ed from galaxy counts. H.E.S.S. observations 
of the blazars H2356-309 and lES 1101-232 bv lAharonian et all J2006h also suggest that the 
CIR B can not be much high er than previously determined lower limits from galaxy counts. 
The lAharonian et al.l (120061 ) Pl.O Extragalactic Background Light (EBL) model, used for 
estimating 7-ray attenuation gives 11, 18 and 16 kJy sr~^ at 1.25, 2.2 and 3.5 fim (26, 25 
and 14 nW m~^ sr~^) and is shown as the top curve in Figure Assuming a power law 
blazar spectrum {dN/dE [photons cm~^ s~^ TeV~^] oc E^^) this model results in a power 
law index for the spectrum of the blazar lES 1101-232 of F = -0.1. The model must be 
scaled down by a factor of 0.45 to the P0.45 model (bottom curve F i gure [T^ to give a power 
law index of at least 1.5, which is considered by lAharonian et al.l (|2006l ) to be the lowest 
acceptable value. The middle curve in the Figure is P0.64 which is t he Pl.O model scaled 



to fit the 10.8ti:} kJy sr'^ {Q-Otll nW m'^ sr'^) val ue reported her e. MapeUi et all (120061 1 



show however that an EBL model, also based on the iWrightl (l2001al ) values at 1.25 and 2.2 
/im, but with a steeper decline from 4 to 10 /xm results in a power law index of F = +0.5. 
They consider F = 0.6 to be the lowest acceptable index based on physical considerations 
and suggest then that while the DIRBE minus 2MASS CIRB values at 1.25 and 2.2 /^m 
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from IWrightl (l2001a] ) and hence the 3.5 fim values reported here and elsewhere in Table [5] do 
require a hard spectrum, and the lower limits from galaxy counts are favored, they can not 
be ruled out based on the current H.E.S.S. data. 



VassilievI (120001 ) showed that it is difficult to untangle attenuation due to the EBL 
from unknown intrinsic properties of blazars, including the total blazar flux or an intrinsic 
deviation from a pure power law spectrum. He showed also that the increased sensitivity of 
H.E.S.S., VERITAS and MAGIC, leading to detections of more blazars at more redshifts now 
and in the near future, and the expanded energy range of these instruments compared with 
their predecessors are crucial to efforts to decouple those effects an d extract information 



about both the EBL and TeV blazars. The MAGIC collaboration (iTeshima et al.l 120071 ) 
claims a detection of 3C 279 at a redshift of 0.5. If we are able to detect TeV sources at such 
large distances, we may be limited to very few CIRB photons as the interaction between 7-ray 
and CIRB photons can put non-trivial limits on the spectrum of the EBL. However, given the 
convergence toward the higher CIRB intensities measured through galaxy counting with the 
faint fringes of the galaxies included and the direct measurements in Table |5l the intensity 
of the CIRB at 3.6/im may be seen as an important anchor for EBL models which will 
allow exploitation of this interaction in determining intrinsic blazar properties and perhaps 
probing certain aspects of fundamental physics such as a boost to the propagation of TeV 
photons by axion production. While the implications of 7-ray attenuation by CIRB photons 
are still limited by the small number of observed sources, this indirect, independent probe of 
the CIRB will improve as more blazars are observed by H.E.S.S., VERITAS and MAGIC. 

Taking this new value from profile-fitting of 10.8^^ ;[ kJy sr^^ (9.0^'^og nW m^^ sr^^) to 
be the fiducial value for the contribution of galaxies to the CIRB, we are left with only a 1.2a 
difference between the values determined via galaxy counts and those determined via the 
direct detection method of "DIRBE minus 2MASS" and others in determining the L-band 
CIRB. Even the more conservative aperture counts are now only 1.8a lower than the direct 
measurements. These numbers are signifi cantly less that the 2.8a discrepancy between direct 



measurements and the lFazio et al.l (120041 ) lower limit of 6.5 kJy sr~^ (5.4 nW m~^ sr~^). It is 
also important to note, however, that the results from profile-fitting do not strongly constrain 
the intensity of the CIRB on the bright end. The probability distribution in Figure M is non- 
Gaussian and the width where the probability drops from its maximum by a factor 
already allows a CIRB of 17.1 kJy sr^^ (14.2 nW m~^ sr^^). Galaxy counts using both 
aperture photometry and profile-fit photometry still bound the CIRB from below, but even 
with a restricted broken power law model for the counts, we are able to find solutions where 
the CIRB goes to large intensities and it is unnecessary to invoke an additional component 
with large numbers at unobservably faint fiuxes. Gamma-ray attenuation may still favor a 
lower intensity, but we do seem to be moving toward a consistent picture of the 3.5-3.6 fim 
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CIRB from direct detection and galaxy counting. The goal then will be to continue to drive 
down the uncertainties in both the galaxy counting and direct measurements. More data 
from Spitzer can help here. One of the limitations of this study is that we have wide shallow 
surveys, but the deep images only survey a tiny solid angle. Multiple deep wide surveys would 
allow more robust galaxy counts at the important magnitudes 16 < m < 20. At 3.6 /im 
these can be done during the post-cryogen extended Spitzer mission when competition with 
other instruments for telescope time will be eliminated. This method is limited, however, 
by confusion not only due to the IRAC pixel size, but more fundamentally due to the real 
confusion between the faint fringes of galaxies in the deep sky. The uncertainties inherent 
in all of these methods means it will likely take an improvement of the zodiacal light models 
to make to make a significant improvement on the precisio n of th ese n ieasurements. The 4 



kJy sr ^ (3.4 nW m ^ sr ^)difference between the IWrightl (Il998l ) and iKelsall et al.l (Il998l ) 



models indicates the level of uncertainty in the zodiacal light s ubtraction. A space m ission 



to directly observe the zodiacal light, of the type suggested in iLevenson et al.l (120071 ) may 
be the only way reduce the dominant source of uncertainty in measuring the intensity of the 
Cosmic InfraRed Background. 



This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under contract 
with NASA. Support for this work was provided by NASA through an award issued by 
JPL/Caltech. 



REFERENCES 

Adelman-McCarthy, J. et al. 2007, ApJ, 172, 634A 

Aharonian, A. G. et al. 2006, Nature 440 1018-1021 

Bertin, E. & Arnouts, S. 1996, A&AS 117, 393B 

Boggess, N. W. et al. 1992, ApJ, 397, 420 

Binder, K. 1997, RPPh, 60, 487b 

Blanton, M. et al. 2001, AJ, 121, 2358 

Bond, J., Carr, B. & Hogan, C. 1986, ApJ, 306, 428 

Cambresy, L., Reach, W., Beichman, C. & Jarrett, T. 2001, ApJ, 555, 563 



- 16 - 

Cole, S. et al. 2001, MNRAS, 326, 255C 

Dickinson, M., Giavalisco, M. & GOODS Team 2002, iastro-ph/0204213 

Dwek, E. & Arendt, R. 1998, ApJ, 508, L9 

Dwek, E., Krennrich, F. & Arendt, R. 2005, ApJ, 634, 155 

Eisenhardt, R, et al. 2004, ApJS, 154, 48E 

Fazio, G., et al. 2004, ApJS, 154, 39F 

Folkes, S. et al. 1999, MNRAS, 308, 459 

Gorjian, V., Wright, E., & Chary, R. 2000, ApJ, 536, 550 

Kelsall, T., et al. 1998, ApJ, 508, 44 

Kochanek, C. et al. 2001, ApJ, 560, 566K 

Kron, R. G. 1980, ApJ43, 305 

Levenson, L., Wright, E., Johnson, B. 2007, ApJ, 666, 34 
Lin, H. et al. 1996, ApJ, 464, 60 

Mapelli, M., Salvaterra, R. & Ferrara, A. 2006, New Astronomy 11 420-430 

Matsumoto, T. 2000, The Institute of Space and Astronautical Science Report SP No. 14, 
p. 179 

Matsumoto, T. et al. 2005, ApJ, 626, 31 
Teshima, M. et al. 2007, larXiv:0709.1475V l 
Simard, L. et al. 2002 ApJS, 142, IS 
Smith, A., 1991 RSPTA, 337, 369S 
Vassiliev, V. 2000 APh, 12, 217V 
Wright, E. 1998 ApJ496: 1-8 
Wright, E. & Reese, E. 2000, ApJ, 545, 43 
Wright, E. 2001a, ApJ, 553, 538 



Wright, E. 2001b, ApJ, 556, L17 



Wright, E. & Johnson, B. 2001, astro-ph/0107205 v2 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 



-18- 



Table 1. DIRBE Dark Spot Coordinates 



Region 


1 


b 


A 


/3 


RA 


DEC 


1 


92.9 


+70.2 


181.9 


+49.5 


13h45m33.6s 


+43d35m59.8s 


2 


318.4 


+76.9 


189.3 


+19.7 


13h05m49.7s 


+14d25mll.9s 


3 


108.9 


+46.4 


153.0 


+70.7 


14h36m23.2s 


+67d46ml8.0s 



Note. — Galactic and ecliptic (J2000) coordinates are given in degrees. 
Equatorial coordinates are J2000. 
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Table 2. Uncertainties and Completeness 



Central Mag^ 




Completeness 


10.5 


± 141% 


100 ± 5% 


11.0 


± 59% 


100 ± 5% 


11.5 


± 59% 


100 ± 5% 


12.0 


± 48% 


100 ± 5% 


12.5 


± 30% 


100 ± 5% 


13.0 


± 18% 


100 ± 5% 


13.5 


± 10% 


100 ± 5% 


14.0 


± 7.5% 


93 ± 10% 


14.5 


± 9.0% 


95 ± 3% 


15.0 


± 11.0% 


83 ± 12% 


15.5 


± 8.9% 


88 ± 5% 


16.0 


± 4.0% 


88 ± 12% 


16.5 


± 9.5% 


74 ± 10% 


17.0 


± 8.3% 


73 ± 10% 


17.5 


± 6.3% 


68 ± 17% 


18.0 


± 7.0% 


71 ± 25% 


18.5 


± 6.6% 


47 ± 15% 


19.0 


± 6.9% 


43 ± 11% 


19.5 


± 7.4% 


37 ± 18% 



'^Magnitudes shown here are profile-fit mag- 
nitudes. 
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Table 3. Most Probable Parameters 





Profile-Fit 


Aperture 


A 


2.3 X 10-^ 


0.5 X 10-^ 


a 


1.85 


2.06 


a 


1.40 


1.75 




17.2 


16.6 
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Table 4. CIRB Corrections [kJy sr ^] 



Photometry 




'Y^{n{m) / comp{m)) f {m) 


Y.'f^modei{fn)f{m) 


Profile Fit 
Aperture 


9.3 

5.3 


12.9 
7.1 


10.8+^} 
q 2+1-2 



Note. — The non-Gaussian nature of the fiux errors, seen in Figures [HI and [T5| 
makes a determination of the errors in the first two numerical columns difficult, 
making necessary the MCMC simulation that provides the above la confidence 
limits on the modeled CIRB. 

'^For comparison, ul^, [nW m"^ sr^^] = jr^Iu [kJy sr^-*^] 



-22- 



Table 5. Previous determinations of the L-band*^ CIRB [kJy sr 



-lib 



Authors 



CIRB 



Zodi-ModeF 



Dwek fc Arendt (1998) 
Goriian et al. f2000l 
Wright k Reese f2000l 
Wright k Johnson (2001) 
Matsumoto (2005) 
Levenson et ah (2007) 
Fazio et al (2004) 
This Work 



11.6 ± 3.4 

12.8 ± 3.8 
14.4 ±3.7 
16.1 ± 4 

16.9 ± 3.5 
15.6 ± 3.3 

> 6.5 

10.8+2;} 



Kelsall et ah (1998) 
Wright (1998) 
Wright (1998) 
Wright (1998) 

Kelsall et al (1998) 
Wright (1998) 

N.A. (galaxy counts) 

N.A. (galaxy counts) 



^Foreground subtraction intensities were measured in the 
DIRBE 3.5 /im band while galaxy count intensities were mea- 
sured in the Spitzer 3.6 /xm band. These results do not distin- 
guish between these nearly identical bands. 



Tor comparison, i>Ii, [nW m ^ sr 



A[/tm 



lu [kJy sr 



-11 



^The iKelsall et aP Jl998h zodiacal light model gives a CIRB 
4.0 kJy sr~^ higher in the L band. 
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Fig. 1. — Galaxy Counts; Aperture vs. Profile-Fit: Galaxy counts vs. magnitude using both 
SExtractor's mag_auto aperture photometry (dashed) and GIM2D profile fit photometry 
(solid). Profile fit magnitudes are approximately 0.5 mag brighter on average than aperture 
magnitudes for the same sources. 
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Fig. 2. — Galaxy Intensity; Aperture vs. Profile-Fit: The Intensity contribution to the CIRB 
from the above galaxy counts is ~ 1.6 times brighter when profile fit magnitudes are used. 
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Fig. 3. — Profile fit counts vs. magnitude in regions 1 2 & 3 (Rl, R2, R3) where beta is the 
echptic latitude and b is the galactic latitude of the region center. Fractional variance in the 
total counts per steradian between the three regions is less than 1%. The fact that the total 
galaxy counts do not show a dependence on galactic latitude, b, indicates successful removal 
of stellar sources in these three regions and supports the above assumption of a csc|b| scaling 
for the star counts. 
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Deep Galaxy Counts 



GOODS-North 



Groth Strip 




16 18 
Vega-Rel Magnitude 



20 



Fig. 4. — Averaged profile-fit galaxy counts from the three regions along with deep counts 
from the EGS and GOODS North fields using the same source detection and photometry. 
Boldness indicates the magnitude range over which each data set was used in compiling the 
total observed galaxy counts. 
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1 



Fig. 5. — The left hand frame shows two of the 40 copies of galaxy A that were inserted into 
the science image for completeness testing. Galaxy B was inserted at the same 40 randomly 
generated coordinates. Both galaxies have profile-fit magnitudes of m = 16.5. However, 
galaxy A has a peak surface brightness of 207.0 kJy sr~^ and is only detected 60% of the 
time while galaxy B has a peak surface brightness of 259.2 kJy sr~^ and is detected 88% of 
the time. 
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Fig. 6. — Completeness vs. peak surface brightness for galaxies at m = 16.5. 
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Fig. 7. — The inserted copies of galaxy A in the left hand frame are the same as those in the 
left panel of Figure [51 Here, the right hand frame shows galaxy C which has a similar peak 
surface brightness to galaxy A (211.9 and 207.0 kJy sr~^ respectively) and both still have 
profile-fit magnitudes of m = 16.5. However, galaxy C has a higher ellipticity and smaller 
half-light radius and while galaxy A is detected 60% of the time, galaxy C is detected 78% 
of the time. 
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Fig. 8. — Results of inserting and re-fitting the indicated number of galaxy images. New F 
is the re-fit flux and Actual F is the total of the pixel values actually inserted. Indicated 
are the number of sample galaxies and the number of random positions at which they were 
inserted. Vertical lines represent the edges of the 0.5 magnitude bins. 
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15 
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Fig. 9. — Results of the MCMC simulation; the distribution of relative probabihties for 
models resulting in a given CIRB. The most probable model from profile-fitting results in a 
CIRB of 10.8^?;} kJy sr-^ (9.0^^;^ nW m ^ sr where the asymmetric range is given by 
the full width where the probability drops from its maximum by a factor e~2. 



-32- 



3.0 




0.5 [ , , , ^ , , , ^ , , , ] 

1.6 1.8 2.0 2.2 

alpha 



Fig. 10. — 68%, 90% and 99% confidence limits in the alpha-beta plane. The dotted line 
shows the convergence prior {a — P < 1) and the cross indicates most probable values. 



-33- 




Fig. 11.— 68%, 90% and 99% confidence limits in the CIRB(m=16.5)-full CIRB plane. The 
CIRB at m = 16.5 is clearly constrained, while the full CIRB is not. The cross indicates the 
most probable values. 



—1 .5 r I I I 

15 16 17 18 19 

Break Magnitude 

Fig. 12. — 68%, 90% and 99% confidence limits on the faint end slope in the power law model 
of the galaxy counts vs. the magnitude at which the break occurs. The cross indicates the 
most probable values. 
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Fig. 13. — Using the averaged counts from our three regions down to m = 16, averaged counts 
from the Groth Strip and GOODS North from 16.5 < m < 17.5 and deep GOODS North 
counts for m > 18, we obtain completeness corrected observed profile-fit galaxy counts at 
each magnitude at 3.5 nm (solid). The shaded area indicates la error bars on the corrected 
counts. The dashed curve shows the modeled "observed" counts. Dot-dashed curve shows 
the most probable model for the "real" counts. 
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Fig. 14. — Using the maximum likelihood model for the "real" galaxy counts (dot-dashed), 
we obtain a total intensity contribution from galaxies as a function of magnitude at 3.5 /im 
(solid). 
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Fig. 15. — Histograms, arranged as in Figure [HI but here, New F are SExtractor aperture 
fluxes and the magnitude at which the histograms are centered are shifted one half-magnitude 
fainter since the galaxies used were selected by their profile-fit magnitudes. 
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Fig. 16. — Scattering the model counts, now with the aperture photometry histograms in 
Figure [151 the MCMC analysis results in a CIRB of 9.2+^;^ kJy sr-^ (7.6lJ;^ nW m'^ sr'^), 
where the asymmetric range is given by the full width where the probability drops from its 
maximum by a factor e~ 2 . The profile-fit results are shown in the background for comparison. 
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Fig. 17. — Curves are the same as in Figure [T3|, but here all counts are determined using 
aperture photometry. 
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Fig. 18. — Profile-Fit vs. Aperture photometry. Top Panel: Average of the uncorrected 
observed counts in regions 1, 2 & 3 as seen in Figure [H Middle Panel: Completeness 
corrected observed counts as seen in Figures [TSlfcfTTl Bottom Panel: Modeled "real" counts 
as determined by the MCMC simulation using both profile-fit and aperture photometry. 
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Fig. 19.— The filled triangle is the 10.8t?i kJy sr ~^ fQ.OtnQ nW m ^ sr ^) value reported 
here. Squares ar e the CIRB values reported in iLevenson et al.l (120071 ). Gr ay stars are 



Matsumotd (120051 ) values from IRTS observations. The open gray triangle is the iFazio et al. 
( I2OO4J ) lower limit at 3.5/x m. For comparison, the upper and lower black curves are the Pl.O 
and P0.45 models used bv lAharonian et al.l (20061) to est i mate the attenuation of TeV 7-rays 
by t he CIRB. Pl.O was norm aliz ed bvlAharonian et al.l (120061 ) to fit the 1.25-3.5 /xm values 
from lDwek &: Arendtl (119981 ) and IWrightl (j2001al ). P0.45 is the Pl.O scaled down by a factor 
of 0.45, which was required, using this shape for the CIRB, to give blazar spectra with power 
law spectral indices of at least 1.5.. Also shown is P0.64 which fits the value reported here. 



